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The dynamical charge fluctuations have been studied in ultra-relativistic heavy-ion collisions by 
using hadronic model simulations, such as UrQMD and HIJING. The evolution of fluctuations has 
been calculated at different time steps during the collision as well as different observation window in 
pseudorapidity (Ar/). The final state effects on the fluctuations have been investigated by varying 
Ar/ and the time steps with the aim of obtaining an optimum observation window for capturing 
maximum fluctuations. It is found that Ar/ between 2.0 and 3.5 gives the best coverage for the 
fluctuations studies. The results of these model calculations for Au+Au collisions at ^/snn = 7.7 
to 200 GeV and for Pb+Pb collisions at 2.76 TeV are presented and compared with available 
experimental data from the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider 
(LHC). 

PACS numbers: 25.75.-q,25.75.Gz,25.75.Nq,12.38.Mh 


I. INTRODUCTION 


The primary goal of heavy-ion collisions at ultra- 
relativistic energies is to explore the signatures of 
the de-confined state of mater, the Quark-Gluon- 
Plasma (QGP). Dedicated experiments at the RHIC 
at Brookhaven National Laboratory (BNL) and the 
LHC at CERN have been setup for studying the 
QGP matter at the high temperatures (T) and low 
baryon chemical potentials ( pb )- Several signatures 
for studying the phase transition from hadronic mat¬ 
ter to QGP have been proposed and as well as stud¬ 
ied in dedicated experiments at BNL and CERN for 
the last few decades. Event-by-event fluctuations of 
conserved charges in limited phase space have been 
widely accepted as one of the most tantalizing sig¬ 
nals of the QGP formation and also for the search of 
the QCD Critical Point EH3. With their large cov¬ 
erage, the STAR experiment at RHIC [s] and the 
ALICE at LHC [9] are ideally suited for the detailed 
study of the QGP matter on an event-by-event ba¬ 
sis. The dynamical charge fluctuations have been 
reported by these experiments [M3. Recent re¬ 
sults from ALICE have shown significant reduction 
in the ratio of charge fluctuations per entropy at the 
LHC energy |l2], confirming to the QGP formation 
in heavy-ion collisions. 

Event-by-event fluctuations of conserved quanti¬ 
ties such as net electric charge and net baryon num¬ 
ber act as distinct signals for the transition from 
hadronic (confined) phase to QGP (de-confined) 
phase. The amount of charge fluctuations is pro¬ 
portional to the squares of the charges present in 


the system, which depend on the state from which 
the charges originate. The system passing through 
a QGP phase has quarks as the charge carriers 
whereas for a hadron gas (HG) the charge carriers 
are the charged hadrons. Thus the charge fluctua¬ 
tions in case of QGP with fractional charges should 
be significantly lower than the HG where the charges 
are integral. Due to the differences in degrees of free¬ 
dom of the two phases, QGP and HG, the magnitude 
of the charge fluctuations are very different. It is es¬ 
timated that for the QGP, charge fluctuations are 
much smaller than the HG [5, [b]. Here the question 
aries whether these primordial fluctuations, either 
from a QGP or from a HG survive during the course 
of the evolution of the system p~6l4l~9l | . The fluctua¬ 
tions observed at the freeze-out depend crucially on 
the equation of state of the system and final state ef¬ 
fects. Non-equilibrium studies at the early partonic 
stage show that large charge fluctuations survive if 
it is accompanied by a large temperature fluctua¬ 
tions at freeze-out [20]. In reality, the measurement 
of charge fluctuations depends on the observation 
window, which is to be properly chosen so that ma¬ 
jority of the fluctuations are captured without being 
affected by the conservation limits p~7l4l9| . 

We have studied the event-by-event dynamical 
net-charge fluctuations originating from the purely 
hadronic state using UrQMD [2l], [22j and HI¬ 
JING [23[ event generators at different times during 
the evolution of hadronic interaction. The dynami¬ 
cal charge fluctuations have been estimated at differ¬ 
ent time steps, and by varying the pseudo-rapidity 
window (Ar/) of the measurement. The main focus 
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is to understand the effect of final state effects which 
diffuse the charge fluctuations at different time and 
A 77 window. Assuming hadronization and freeze- 
out occur roughly at 5 fm/c and 30 fm/c, respec¬ 
tively, we have calculated the fluctuations of the sys¬ 
tem at 5 fm/c, 30 fm/c, and at a much latter time of 
100 fm/c, where all possible interactions must seize. 

This paper is organized as follows. The measure of 
dynamical charge fluctuations in heavy-ion-collisions 
are discussed in Section-II. In Section-Ill, we present 
particle multiplicity distributions at Arj = 1, for dif¬ 
ferent time steps for central Au+Au collisions at 
yj snn = 200 GeV, using UrQMD. The measure of 
dynamical charge fluctuations at different time steps 
and A ?7 window are discussed in Section-IV. The re¬ 
sults of the calculations from hadronic models are 
presented in Section-V, along with the experimental 
data from STAR and ALICE. The paper is summa¬ 
rized in Section-VI. 


II. NET-CHARGE FLUCTUATIONS 


The Net-charge and the total charge of a sys¬ 
tem are denoted in terms of Q = N+ — N _ and 
iVch = N + + AA, where N + and AA are the mul¬ 
tiplicities of positive and negative charged particles, 
respectively. The net-charge fluctuations can be ex¬ 
pressed in terms of its ratio to entropy in order to 
take the volume term into account. Thus, one of the 
observables for net-charge fluctuations is Q: 


D = 


A_m 

N ch ’ 


(1) 


where SQ 2 is the variance of the net-charge. The 
value of D has been estimated by theoretical models 
for a QGP and a HG by taking various final state 
effects into account [Els HE HE HEHE ■ Early es¬ 
timations had put the value of D to be approxi¬ 
mately 4 times smaller for a QGP compared to a 
HG. For a HG, resonance decays including those of 
neutral particles introduce additional correlation be¬ 
tween charged particles, which reduces the value of 
d HE HE. Present understandings put the value of 
D to be 1 — 1.5 for a QGP and 2.8 for a HG. In 
all cases, the signal gets diffused from hadronization 
time to freeze-out because of the final state interac¬ 
tions which need to be taken into account EE HE- 
Net-charge fluctuations, measured in terms of D , 
have contributions from statistical as well as dynam¬ 
ical origin. It is a rather difficult task to estimate the 
dynamical component from the total fluctuations. A 
novel method of estimation of the dynamical fluctu¬ 
ations has been proposed, which takes into account 

the correlation strengths between +-(-,-and H— 

charged particle pairs [IE- The difference between 
the relative number of positive (N + ) and negative 


(AA) charged particles can be expressed in terms of 
its second moment as, 


A— = 




( 2 ) 


Here, the notation “( }”denotes average over the 
ensemble of events. Assuming independent parti¬ 
cle production mechanism, the value of v _| in the 

Poisson limit can be expressed as, 

l/+_ ’ BtBt = (A+) + (AA)‘ ® 

The dynamical component is then evaluated as a 
difference between the two measured fluctuations, 
expressed as, 


A-— ,dyn — A— A —,stat- (4) 


This can be expanded as, 


A—,dyn 


(N + (N + - 1)} 
(N + y 


(AA(AA — !)) 
<■ N -) 2 
(AAAA) 

~ 2 (n + )(n_)- 


( 5 ) 


A stronger correlation between H— pairs compared 
to +- 1 - and-pairs yields a negative value of 

A—,(Un¬ 
it can be seen that the a —,dyn is related to the 
net-charge fluctuations, D by, 


(-Nch)A—,dyn = D- 4. (6) 

By determining _ jC j yn in the experiments, one can 

have access to net-charge fluctuations. 

The magnitude of net charge fluctuations is lim¬ 
ited by global conservation of charged particles HE ■ 
Considering the effect of global charge conservation, 
the dynamical fluctuations need to be corrected by 
a factor of A--,d yn = —i/{N^), where (N^) is the 
average of total number of charged particles pro¬ 
duced over full phase space. The corrected value 
of a _,dyn after considering the global charge con¬ 

servation and finite acceptance is 

*4-4 y„ = ^+-,dyn+^. (7) 

iV4 7r 

The modified value of the net-charge fluctuations 
turns out to be: 


D = (A ch )i/+°" dyn + 4. ( 8 ) 

In the rest of the article we will evaluate 
(N c ^-Ayn and D for different center-of-mass en¬ 
ergies. 
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III. MULTIPLICITY DISTRIBUTIONS AT 
DIFFERENT TIME STEPS 

In order to understand the evolution of multi¬ 
plicity distributions of different particle species at 
different time steps, we have used UrQMD model 
simulations for Au+Au collisions corresponding to 
RHIC energies. The UrQMD model simulates the 
microscopic transport of covariant propagation of 
quarks and di-quarks with hadronic degrees of free¬ 
dom. The formation of hadrons is introduced by the 
color string fragmentation. Various resonances and 
their decay along with re-scattering among hadrons 
have been incorporated during the evolution [22^ . 
This model helps to explore the evolution of con¬ 
served charge fluctuations and their distribution at 
different time steps in the hadronic medium. 

In the present study, the UrQMD model has 
been used to simulate Au+Au collisions at vari¬ 
ous collision energies. The event-by-event distribu¬ 
tions of different charged particle and anti-particle 
species are estimated at time 5 fm/c, 30 fm/c and 
100 fm/c after the collision. Multiplicity distribu¬ 
tions within 1171 < 1.0 and transverse momentum 
range of 0.2 < pp < 5.0 GeV/c are presented in 
Fig. [T| for central (0 — 5% centrality) Au+Au colli¬ 
sions at ^/snn =200 GeV Multiplicity distributions 
of charged particles (N + and 1V_), pions (7r + and 
7T“), kaons ( K + and K~), and protons (p and p) 
are shown for the three time steps. The distributions 
shift to the right as the system evolve with time in 
going from 5 fm/c to 30 fm/c and 100 fm/c. The 
shifts for the N+ and N -, pions and kaons are quite 
appreciable, whereas protons and anti-protons are 
less affected. The multiplicity distributions of pions 
and kaons mainly contribute the change in total pos¬ 
itive and negative charged multiplicity distribution. 
The shift of the kaon multiplicity distributions af¬ 
ter 5 fm/c could occur because the kaon production 
from meson-meson and baryon-meson interactions, 
as implemented in UrQMD model, dominate during 
this time. Additional change at higher multiplicity 
may be due to rescattering and resonance decays in 
a given phase space. Because of their higher masses, 
the distributions for protons and anti-protons com¬ 
pared to those of the pions and kaons, are less af¬ 
fected during the evolution of the system. The pro¬ 
ton number is expected to diffuse more slowly be¬ 
cause of re-scattering Q. This change of multiplicity 
distributions is maximum in larger pseudo-rapidity 
window. 

Due to the final sate effects, the change of the 
shape of multiplicity distributions may affect various 
event-by-event observables. The fluctuations of mul¬ 
tiplicity distributions diffuse at different time scales 
in heavy-ion-collisions, in the rapidity space. Hence, 
it is expected that different fluctuations measures 






FIG. 1: (Color online) Multiplicity distributions for 
Au+Au collisions at U s nn =200 GeV within 1 77 I < 0.5 
and 0.2 < pr < 5.0 GeV/c at different time steps 5 
fm/c, 30 fm/c, and 100 fm/c for (a) positive and nega¬ 
tive charged particles, (b) 7r + and 7r _ , (c) K + and K ~, 
and (d) p and p 
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FIG. 2: (Color online) The values of (N c h)i/_|_,dyn (up¬ 

per panel) and (A r c h)t'“"dyn (lower panel), plotted as 
functions of Ar/ window using UrQMD model at two 
different time steps for central (0-5%) Au+Au collisions 
at ,/^NN — 200 GeV. 


may be affected differently with the time evolution 
in a given phase space. In the next sections, we 
present the dynamical charge fluctuations measures 
at different time steps for Au+Au collisions at 200 
GeV using UrQMD model. 


IV. FLUCTUATIONS AS A FUNCTION OF 

At] 


The evolutions of t'+°V dyn are studied by us¬ 
ing UrQMD by varying different A?y windows for 
different time steps. The main goal of this exer¬ 
cise is to understand the evolution of fluctuations 
through purely hadronic medium as well as to find 
an optimum coverage where most of the fluctuations 
can be measured. This information helps to under¬ 
stand the evolution of fluctuations through purely 
hadronic medium, as charge fluctuations are sup¬ 
posed to be diffused with the increase in Ar) win¬ 
dow. Total charge of a system is conserved leading 
to vanishing net-charge fluctuations for full cover¬ 
age. At the same time studying fluctuations in a 
very small At] window may not be ideal for cap¬ 
turing most of the initial fluctuations. An optimum 
coverage is to be obtained by taking these into ac¬ 
count. 

To obtain the optimum value of fluctuations, tak¬ 
ing all effects into account, we have considered 
A T] range from 0.2 to 10.0. The fluctuations are 
calculated for central (0-5%) collisions. To avoid 
the dependence on the central bin width, the value 

of iq_dyn is determined using unit bin method. In 

this method, value of z+_ ^dynim) for each multiplic- 



FIG. 3: (Color online) The ratios of {N c h) +_ ,dyn (up¬ 

per panel) and (lV c h) Vj°.I( dyn (lower panel), with respect 
to their are normalized values at the smallest Ar/ of 0.2 
for central Au+Au collisions at ^/snn=200 GeV. The ra¬ 
tios are plotted as a function of A^for two different time 
steps, 5 fm/c and 30 fm/c. 


ity is calculated and then averaged over the width of 
particular centrality with the weights corresponding 
to relative cross section. The weighted average for 
i'-l_ t dyn are calculated as: 


V-\ — ,dyn\JYlmin 


TCI <C Mmax ) — 


E v+-,dyn{m)p(m ) 


E p( m ) 


(9) 


Here, p(m) is the weight of particular centrality m. 

Finally, the corrected values of z+_ jdyn have been 

obtained using Eq. (7). 

Figure [2] shows both the uncorrected 
((iV ch ) Z z + _, dyn ) and corrected ((V ch )z/™" dyn ) 
values of fluctuations as a function of Ar] win¬ 
dow for central (0-5%) Au+Au collisions at 
■y/SNN = 200GeV, obtained from UrQMD. The 
results are presented for two time steps, 5 fm/c, 
and 30 fm/c. The trends for the uncorrected and 
corrected values of fluctuations are observed to 
be very different. The upper panel of the figure 

shows that _ ,dyn keep decreasing with 

the increase in Ar], This is unphysical as the 
fluctuations should vanish for measurements at 
the full coverage. The nature of (N c h)^™I r dyn , 
on the other hand, shows different trend, where 
the values decrease up to Ar] values of 2 to 2.5, 
then remain constant till about Ar] =3.5, and then 
increase as per the expectations. The values of 
<iv ch K°v dyn tend to zero at the highest Ar] due 
to the global charge conservation. This decreasing 
trend of (IV c h)i'+ 0 ” dyn up to Ar]~2 is due to the 









5 


strengthening of multiplicity correlations with 
increase in A rj. 

The nature of the fluctuations, at two time steps, 
as a function of A ry may be better understood by 
plotting the ratio of the fluctuations at different 
Ary values with respect to a particular Ary (normal¬ 
izing with respect to smallest Arj). Figure [3] shows 
the ratios of (N ch )v+-,d y n and (V ch )r/j_°V dyn with 
respect to their values at Ar) = 0.2. From the 
upper panel of the figure, it is seen that the un¬ 
corrected normalized (N c h)z+_^yn ratios increase 

monotonously with increase in Arj. On the other 
hand, the corrected normalized N c hV < ^L r dyn values 
increase up to Ary~2, then remain constant up to 
Ar] =3.5. As the hadronic system evolves, it en¬ 
counters more and more rescattering and resonance 
decay as compared to smaller pseudo-rapidity win¬ 
dow. Within Ar] range of 2.0 and 3.5, the diffusion 
of dynamical charge fluctuations may remain insen¬ 
sitive. Going beyond Ar] = 3.5, the fluctuations de¬ 
crease due to the dilution of correlations and effect 
of global charge conservation. Near Ary of 8.0, the 
fluctuations are close to zero. Going to higher Ary, 

the ratio goes below zero, indicating v d _ !( j yn becomes 

positive. This could happen because -l-ve and —ve 
charged particles become uncorrelated, possibly be¬ 
cause of inclusion of spectator particles. In addition 
to the dependence of fluctuations on Ary, Fig. [3]also 
gives the time dependence of fluctuations for wide 
Ar] windows compared to a narrow bin of Ary = 0.2. 
It is observed that the fluctuations for wide Ary com¬ 
pared to the corresponding narrow Ary are more 
pronounced at a time of 5 fm/c compared to the 
corresponding values at a later time of 30 fm/c. 

From the present study, we conclude that the op¬ 
timal coverages for observing the charge fluctuations 
are for Ary =2 — 3.5 for ^/snn = 200 GeV. For lower 
energies, the Ary window will be somewhat lower. 
These values are in confirmation with earlier pub¬ 
lished results EMI- The Ary dependance of charge 
fluctuations may give information about the proper¬ 
ties of the hot and dense medium created in heavy- 
ion collisions 0- 


V. COMPARISON OF MODEL 
CALCULATIONS WITH EXPERIMENTAL 
DATA 


Net-charge fluctuations have been measured 
by experiments at CERN-SPS, RHIC and LHC. 
Recently ALICE experiment published the net- 
charge fluctuations for Pb+Pb collisions at 
yj snn =2.76 TeV 0- The results from the STAR 
experiment at RHIC energies had been published 
earlier 0. The measured values of net-charge 
fluctuations are presented in Fig. [I] where both 
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FIG. 4: (Color online) (A 0 h)i'+ 0 T ( jyn (left-axis) and cor¬ 
responding values of D (right-axis) as a function center 
of mass energy in Au+Au or Pb+Pb collisions from HI- 
JING and UrQMD event generators for different Ar] win¬ 
dows. Estimations for fluctuations originating from pion 
gas, hadron resonance gas and QGP are indicated. 
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HIJING An = 1.0 

- UrQMD Ail = 1.0 

- UrQMD Ail = 2.0 

- UrQMD Ari = 3.0 

- UrQMD Ari = 4.0 
STAR Au+Au Ari = 1.0 
ALICE Pb-Pb Ar\ = 1.0 
ALICE Pb-Pb At) = 1.6 

Pion Gas 

Hadron Resonance Gas 
QGP 


(N ch )^ dyn and D are plotted as a function of 
center of mass energy for Pb+Pb collisions at LHC 
and Au+Au collisions at RHIC. The STAR results 
are measured for Ar] =1.0 and the ALICE results 
are shown for both Ary = 1.0 and 1.6. The values of 

dynamical net-charge fluctuations, tq_u yn , remain 

negative at all cases that implies the existence of 
finite correlation between +ve and —ve particles. 
The fluctuations are observed to decrease as the 
center of mass energy increases. 

We have calculated the net-charge fluctuations 
for Au+Au collisions at ^/snn = 7.7, 11.7, 19.0, 
27, 39, 62, and 200 GeV, using two hadronic mod¬ 
els, HIJING and UrQMD. The HIJING model is a 
perturbative-QCD inspired model which contains jet 
and mini-jet formation mechanism. On the other 
hand, UrQMD is a transport model which contains 
various resonance decays and, elastic and inelas¬ 
tic interactions. The results are superimposed in 
Fig. |T] The HIJING calculations are performed at 
Ar] =1.0 for 0-5% central collisions. The UrQMD 
results are for time at 30 fm/c, and for a set of val¬ 
ues at Ary from 1.0 to 4.0. The values of D from 
these model calculations are within the pion gas and 
hadron resonance gas limits. 

The net-charge fluctuations obtained from experi¬ 
mental measurements at RHIC energies for Ary = 1.0 
are within the pion gas and hadron resonance gas 
(HRG) limits. Extending the Ary range will be 
advantageous for better understanding the fluctu¬ 
ations. At ^/snn =2.76 TeV corresponding to LHC 
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energy, the result for A7? =1.0 for central collision is 
below the HRG limit. At A rj = 1.6, the fluctuations 
further decrease. The values of D being within the 
HRG limit and QGP imply that at LHC energy the 
fluctuations have their origin in the QGP phase. 

VI. SUMMARY 

We have studied the dynamical charge fluctua¬ 
tions at different time steps using UrQMD model 
for Au+Au collisions. The positive and nega¬ 
tive charged particle multiplicity distributions, at 
A7? =1.0 for central collisions, change with time. 
It is found that contributions at different time steps 
for protons and anti-protons are less as compared to 
those of the pions and kaons. Dynamical fluctua¬ 
tions are studied using i/j°" dyn , corrected for global 
charge fluctuations. The net-charge fluctuations, ex¬ 
pressed in terms of D and (N c h)i'™L T /iyn are stud¬ 
ied for a range of At?, from a narrow window of 
0.2 to the maximum of 10.0 for Au+Au collisions 
at -y/SNN = 200 GeV. One of the major goals of the 
present study is to find an optimum A 7? window 
for which maximum amount of charge fluctuations, 
originating from the early stages of the collision, can 


be captured. We find that with increasing A 7? win¬ 
dow, the value of fluctuations increase, indicating 
final state effects, such as resonance decay and re¬ 
scattering. The value of D does not grow any more 
beyond A77 = 2.0. On the other hand, D remains 
constant till At? = 3.5, and then decreases close to 
zero for A7? = 10.0. This observation confirms the 
charge conservation scenario. From this study, we 
can conclude that the optimum value of charge fluc¬ 
tuations are captured for At? = 2.0 — 3.5. 

The charge fluctuations, obtained from HIJING 
and UrQMD models are compared to the experimen¬ 
tal data at RHIC and LHC energies. It is observed 
that a value of At? around 2.0 is ideal at all energies 
for studying charge fluctuations. As expected, the 
results from the model calculations remain within 
the limit of pion gas and hadron resonance gas val¬ 
ues for all energies. 
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